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Abstract 

Complete description of the singular sectors of the 1-layer Benney 
system ( classical long wave equation) and dToda system is presented. 
Associated Euler-Poisson-Darboux equations E(l/2,l/2) and E(-l/2,-l/2) 
are the main tool in the analysis. A complete list of solutions of the 
1-layer Benney system depending on two parameters and belonging to 
the singular sector is given. Relation between Euler-Poisson-Darboux 
equations E(e, e) with opposite sign of e is discussed. 



1 Introduction 

The 1-layer Benney system (classical long wave equation) 

u t + uu x + v x = 0, 

(1) 

v t + (uv) x = 
and dToda equation v xx = (logv)tt or equivalcntly the system 

u t + v x = 0, 

(2) 

^ v t + vu x = 

are the two distinguished integrable systems of hydrodynamical type (see e.g. 
[1, 2]). 1-layer Benney system describes long waves in shallow water with free 
surface in gravitational field. It is the dispersionless limit of the nonlinear 
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Schrodinger equation [3]. Recently, the 1-layer Benney (.6(1)) system became 
a crucial ingredient in the analysis of the universality of critical behavior for 
nonlinear equations [4]. The dToda equation is the 1+1-dimcnsional version of 
the Boyer-Finley equation from the general relativity [5] . It shows up in various 
problems of fluid mechanics ( see e.g. [6]). It is known also (see e.g. [7]) that the 
hodograph equations of the dToda hierarchy determine the large TV-limit of the 
Hcrmitian model in random matrix theory. In general, these two systems are an 
excellent laboratory for studying properties of integrable hydrodynamical type 
systems. 

In the present work we analyze the structures of the set of hodograph equa- 
tions of the -B(l) hierarchy and dToda hierarchy in terms of its Riemann invari- 
ants. These hodograph solutions describe the critical points 

dW 

M = > i = 1 ' 2 ' (3) 

of a function W = W(t, fj\, fa) which depend linearly on the coordinates t, 
where t denotes the flow parameters of the Benney hierarchy and dToda hier- 
archy respectively and obey an Euler-Poisson-Darboux equations E(e,e) [8] 

where for the Benney system one has e = 1/2 and for the dToda system e = 
— 1/2. The equation (4) and its multidimensional version are well known for 
a long time in classical geometry [8] . Its relevance to the theory of Whitham 
equations has been demonstrated recently in the papers [9] -[11]. Here we will 
use classical notation E(e, fi) for the Euler-Poisson-Darboux equation proposed 
in [8] where such equations with different e and fi have been studied too. 

If we denote by M. the set of solutions (i,/3) (/3i ^ fa) of the hodograph 
equations (3), we may distinguish a regular and a singular sector in M. 

M=M rc SuM sin g, 

such that given (t, (3) g M. 

('•« < * - "*» (^ww) - " 

The elements of A4 ie &, correspond to the case when the system (3) is uniquely 
solvable and hence, it defines a unique solution (3(t). The singular class M sm & 
represents degenerate critical points of the function W and are the points on 
which the implicit solutions (3(t) of the hodograph equations exhibit "gradient 
catastrophe" behaviour. As we will see in this paper, the Euler-Poisson-Darboux 
equation is of great help to analyze the structure of A4 sm §. As the illustration 
of the general result a complete list of solutions of the 1-layer Benney hierarchy 
from 7VJ sm § depending on two parameters is presented. We also discuss the 
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relation between Euler-Poisson-Darboux equations with opposite a and Euler- 
Poisson-Darboux equations for symmetries and densities of integrals of motion 
for integrable hydrodynamical type systems. 



2 1-layer Benney hierarchy and its singular sec- 
tor 

The B(l) system (1) is a member of a dispersionless integrable hierarchy of 
deformations of the curve (sec e.g. [14, 15]). 

p 2 = {\-p 1 ){\-p 2 ). (5) 

where u — — ((3\ + /3 2 ), v = — P2) 2 ■ The flows (3{t) are characterized by the 
following condition: There exists a family of functions S(\, t, (3) satisfying 

d tn S(X, t, 0(t)) = ft„(A, 0(t)), n > 1. (6) 

where 

W = ( ,,ul-J , >A>-A><*-A). (7) 



where ® denotes the standard projection on the positive powers of A. Functions S 
which satisfy (6) are referred to as action functions in the theory of dispersionless 
integrable systems (see e.g. [16]). Notice that for n = 1 (6) reads 

dS 

x:=h, 

ox 

so that the sytem (6) is equivalent to 

d tn P = d x fi„, (8) 

and, in terms of Riemann invariants /3, it can be rewritten in the hydrodynamical 
form 

/ A™ \ 

The i 2 -flow of this hierarchy is the B(l) system (1) (t := t 2 ) 

(10) 



dtfa = ^(3 02 + 0l)ft x . 

For v > the B(l) system is hyperbolic while for v < it is elliptic. 

It was proved in [12] that the system (3) for the critical points of the function 

W {m -.= - V ^ t] , (11) 

J 7 2 ? 7r v / (A-/3i)(A-/3 2 )' 
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where V(X, t) — J2n>i ^« 1S a system of hodograph equations for the Benney 
hierarchy. Moreover, the action function for the corresponding solutions is given 

by 

S(X,t,0) = Yl tnftn&P) = KX,t,0) v/(A-/3i)(A-/3 2 ). (12) 

n>l 

where 

V(A-ft)(A-ft)J© 

Obviously, the function VF satisfies the Euler-Poisson-Darboux equation £7(1/2, 1/2). 
Written explicitly, W represents itself the series 

W=^(p 1 + ft) + | (3ft 2 + 2ftft + 3ft 2 ) + ^| (5ft 3 + 3ft ft + 3ft ft 2 + 5/3 3 ^ 



2^ 8 v ^ rz ' 16 

*4 /or Q 4 i oo«3o i 1Q o2o2 i onfl fl3 , oro4\ 



2/ 



+ ^ (35ft + 20ftft + 18ft ft 2 + 20ft ft + 35ft 4 ) + • • • . (13) 
The hodograph equations (3) with t n = for n > 5 take the form 
8a; + 4t 2 (3ft + ft) + 3t 3 (5ft 2 + 2ftft + ft 2 ) + ^(UOft 3 + 60ft 2 ft + 18ftft 2 + 20ft?) = 0, 



8a; + 4£ 2 (ft + 3ft) + 3t 3 (ft 2 + 2ft ft + 5ft 2 ) + ^ (140ft? + 60ft 2 ft + 18ftft 2 + 20ft 3 ) = 0. 

(14) 

Detailed analysis of equations (14) will be performed in section 3. Here, we 
would like to make two observations. First, one is that the formulae (14) point 
out on the possible alternative interpretation of the times i 2 , t 3 , t^,... of the 
B(l) hierarchy. Namely, taking i 2 = in the formulae (14), we see that £3 and 
£4 are parameters appearing in the initial data ft(a;,i 2 = 0) and ft(a;,t 2 = 0). 
Thus, one can view hodograph equations (3) (in particular, equations (14)) as 
equations describing the time evolution of the family of initial data for the B(\) 
system , parametrized by the variables £3, £4, t§,... 

Second observation concerns with the elliptic version of the B(l) system. In 
this case ft = ft and the system (10) reduces to the single equation 

d t p=\{3p + ft ft, t:=t 2 , /3:=ft. (15) 
This equation is equivalent to the nonlinear Beltrami equation 

2i - 3/3 - p 



2i + 3/3 + p 



ft, (16) 



where z = x + it. This fact indicates that the theory of quasi-conformal 
mappings (see e.g. [17]) can be relevant for the analysis of properties of the 
elliptic B(l) system (1) (v < 0). Hence, since the elliptic B(l) system is the 
quasiclassical limit [3] of the focusing nonlinear Schrodinger (NLS) equation 

ieipt + y ipxx + M tp = 0, 
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with ib = A exp (- s), u = — (— - ,v = - \ib\ 2 , e -> 0, the quasicon- 
\e / 2i\ ip ib ) 

formal mapping can be useful also in the study of the small dispersion limit of 

the focusing NLS equation (compare with [4]). 

In order to analyse singular sector of the 1-layer Benney hierarchy we first 

observe that due to the Euler-Poisson-Darboux equation for given (t,/3) € A4, 

as a consequence of (4) one has 



d 2 W 

0. (17) 



Consequently 



det 



d fad fa 
d 2 W \ d 2 W d 2 W 



\ _ cr w cr w 

J = ~dW'~dW ( ' 



.dfadfaJ dp 2 dp. 



Thus, we have 

Proposition 1. Given (t,/3) e M then 

1. (t,/3) e M re 9 if and only if ^ ^ Oand^ ^ 0. 

2. (t,/3) G M, s * n 9 if and only at least one of the derivatives ^pr, 
vanishes. 

Furthermore, using (4) it follows easily that at any point (t,/3) e M. all 
mixed derivatives d 1 ^ d 3 ^ W can be expressed in terms of linear combination of 

derivatives d^ ± W and d™ 2 W. Hence if we define M^§ 2 as the set of points 
(t, (3) e M such that 

d n '+ 2 W d k W 
it follows that 

-M sing = IJ M^l, 

ni+n 2 >l 

where 

We may characterize the classes M ni ,n 2 of the singular sector in terms of 
the behaviour of S(X) at A = fa (i = 1, 2) . Indeed the derivative d^ 1-1 W with 
fc > 1 is proportional to the integral 



dA V(X,t) I dA h(X) 



2it: (X-fa) k +^(X-fa)(X-fa) J 7 2i7r(A-ft) 
and a similar result follows for the derivatives d k , +1 W with k > 2. As a conse- 

P2 — 

quence we have 



A=/3i 
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Proposition 2. A point (i,/3) € belongs to the singularity class Mnun 2 if 
and only if 



S(\,t,P) ~ (A- fr)-*- as A ->■ ft, (i = l,2) 



(20) 



3 Explicit determination of singular sectors 

It is easy to see that the singular classes A4^n 2 can be determined by means of 
a system of m + ri2 constraints for the coordinates t. Indeed, the points (t, [3) 

of A4^§ 2 are characterized by the equations 
d k W 

— - F =o, \Jl<k<m + l, i = 1,2, (21) 



9ft 



and 



X5-7^0, i = 1,2. 



(22) 



Now the observation is that the jacobian matrix of the the system of two equa- 
tions 

Qm + l w 

= 0, i = 1, 2 



is not singular as 




(23) 



* 0. 



(24) 



Indeed, we notice that as a consequence of (4) the derivatives outside the diag- 
onal of A are linear combinations of the derivatives {dp. W, 1 < k < rii + 1, i = 
1,2}, so that from (21)-(22) we have 

Qn 1+ 2 W Qn 2 +2 W 

* ^+r^o. 



ni+2 



dp. 



Therefore, one can solve (23) and get a solution (3{t). Substituting this solution 
in the remaining equations (21) gives m + n-i constraints of the form 

/fc(t)=0, k = 1, ...,m +n 2 . 

It is not difficult to determine the solutions of (21)-(22) in two simple cases: 
with one parameter ts {ti = t$ = ■ ■ ■ = 0), and with two parameters t%, ti 
(£5 = te = ■ ■ ■ = 0). We have that in this case 



A^sing = ^smg {J ^sing 



G 



with M^ 1 ^ defined by 

-45£ 4 t| + 180* 2 *i*3 + ^(8*2*4 - 3*i)V*! (3*§-8* 2 * 4 ) 



1. X = 



360*1 

5*3*4 + VlE a/*! (3*1 - 8* 2 * 4 ) fl _ -3*3*4 + \/l5-\AI ( 3i 3 - 8i 2*4) 



20*1 ' P " 12*2 



2. a; 



-45*4*| + 180*2*1*3 - V15(8*2*4 - 3*|) (3*| - 8* 2 * 4 ) 
360*1 

-5*3*4 + \/l5V*l (3*| - 8* 2 * 4 ) „ _ 3*3*4 + VTEy/tl (3*1 - 8*2*4) 



20*1 ' HZ 12*1 

and Mor g b y 



-45* 4 *i + 180*2*1*3 - Vl5(8*2*4 - 3*i)V*| (3*| - 8*2*4) 



360*1 



3*3*4 + \/l5V*4 (3*1 - 8*2*4) „ _ -5*3*4 + \/15a/*4 (3*1 - 8*2*4) 



2 = 



12*1 20 *4 



-45*4*1 + 180*2*1*3 + 715(8*2*4 - 3*|)y/*l (3*1 - 8*2*4) 

360*1 



-3*3*4 + yfiEy/t\ (3*1 - 8*3*4) „ 5*3*4 + 715^1 (3*1 - 8*2*4) 



n = T7T75 1 P2 



12*1 ' 20*1 

4 Singular sector of the elliptic B{1) system 

Now, we will consider the elliptic B(l) system (1). Singular sector M sm & has in 
this case a structure which is quite different from that of the hyperbolic system. 
Indeed, since Pi = P 2 , the function W for real x, *2, *3, . . . is the real valued 
function 

W(t,p,p) = W(t,p,p), 

and the hodograph equation (3) has the form of the Cauchy-Riemann condition 

(P = Pi) 

dW 

^=r = 0. (25) 

Regular and singular sectors M ie & and M sm & are defined as the sets (*, P, p) 
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of solutions of equation (25) such that the hermitian form 



d 2 W 



d 2 W 
~~dj 2 



dp 2 + 2 



d 2 W 



dpdp 



= d[3df3 



d 2 W 



d/3 



is nondegenerate or degenerate, respectively. For unreduced solutions (P ^ P), 
the corresponding Euler-Poisson-Darboux equation implies that 



and, hence 





d 2 W 




dpdp 


d 2 W 


d 2 W 


dp 2 


dpdp 


d 2 W 


d 2 W 


dpdp 


df 



= = 0, 



d 2 W 



op 2 



(26) 



Thus, one has 



Proposition 3. For unreduced solutions of the elliptic B(l) system, the regular 
sector M re 9 is defined by the condition 



dW 

— =■ = 0, 

dp 



d 2 W , 
-=2- ^ 0. 

dp 



(27) 



A similar analysis to that of the hyperbolic case readily leads to 

Proposition 4. Singular sector M sm 9 of the elliptic B(l) system (1) is the 
union of the subspaces Ain"^ {n = 1,2,3,...) defined as 

{ — ■ d k W M n + 2 ) W 1 

M s n mg = (t,P,P) G M sm 9 : = 0, k = l,...,n + l; ^ 



dp 



df n+2) 



(28) 

Solutions belonging to Mm™® are defined on a subspace of codimension 2n in 
the space of parameters x, t 2 , tz, ... . 

So, in the elliptic case, gradient catastrophe happens in the point (x,t) at 

sins? 

fixed parameters t 2 , H,... Similar to the hyperbolic case, the subspace M n 
is not empty if at least n parameters t 2 , t^,...,t n+ i are different from zero in the 
formula (13). 

It is instructive to rewrite the formula (13) for the function W in terms of 
the real and imaginary part of P\. i.e. P\ = U + i V: 

w = xu + t 2 (u 2 - \v 2 ) + t 3 (u 3 - \uv 2 ) + u(u 4 - su 2 v 2 + fv 4 ) 

2 2 8 



+ t 5 (U 5 - 5U 3 V 2 + ^-UV 4 ) + 



(29) 
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This formula explicitly shows the character of elliptic singularities exhibited for 
the function W for various values of parameters t 2 , H, .... 

Basic equations (25), (13) and also conditions (28) defining subspaces M^ n ^ 
can be easily rewritten in terms of the original variables u and v. Since 

dW dW . dW 

+ i *J— v - 



dp du dv 

the hodograph equation (25) becomes (for v ^ 0) 

dW „ dW 

-5- = 0, -5- = 0, 30 

du dv 

while the Euler-Poisson-Darboux equation and equation (30) take the form 

d 2 W d 2 W 



du 2 dv 2 
For the subspace A4^ n ^ conditions (28) are 



= 0. (31) 



Since 



dW n d 2 W n d 3 W , n 

aF = ' (32) 

d 2 W d 2 W . d 2 W d 2 W 1 dW 

2 1 v — v — — — + v ■ 



d /3 2 du 2 dudv dv 2 2 dv 7 

one concludes taking into account equation (30) and (31) that the second con- 
dition (32) is satisfied if and only if 

d 2 W n d 2 W n d 2 W 

= 0, -—- = 0, — — = 0. (33) 



du 2 ' dv 2 dudv 
sins 

Thus, the subspace M. 1 is characterized by the conditions (30), (33) and by 
requirement of nonvanishing third order derivatives of W . 

In order to compare these conditions with those of paper [4] , we first observe 
that the -B(l) system (1) is converted into the system (1.8) by the substitution 
u — y v 1 v — y — u. Then, with the choice 

W = f(u,v) + xv — uvt, 

the hodograph equations (30) become equations (2.4) of [4] and equation (31) 
is reduced to their equation (2.5). Finally, with such a choice, the conditions 
(33) are converted to the condition (2.12) from the paper [4]. 

Finally, we note that according to the proposition 4 for the subspace A4 1 , 
the codimension of the corresponding subspace of (x, t 2 , t 3 , . . . ) is equal to two 
and the function W with t n — 0, n > 4 i.e. 

w = xu + t 2 (u 2 - \v 2 ) + t 3 (u 3 - luv 2 ), 

o Z 

exhibits the elliptic umbilic singularity according to Thorn's classification [20] 
(see also [17]- [19]). These results reproduce those originally obtained in the 
paper [4] (formula (4.2)) 
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5 5. dToda hierarchy. 

Now let us consider the function 



W T (x,0M= J ^V T (x,X)yJ(l-^){l-^) (34) 

where Vr{x, A) = ^ n > X n x n . Critical points for this function are defined by 
the equations 

It is a simple check to see that Wt obeys the Euler-Poisson-Darboux equation 
of the type £(-1/2,-1/2) 

2 ^- fii) mm = ~ i 'W~~dk ) - (36) 

Written explicitly the function Wt is the series 

? 2 )(/3i-/3 2 ) 2 - I ^X3(5^ + 

(37) 



^ = -^o(/3i+/3 2 )-^ 1 (/3 1 -/3 2 ) 2 -^ 



while the hodograph equations take the form 

x + ^i(/3i - ft) + lx 2 (3/3 2 - 2^02 - 0i) + . . . = 0, 
z o 

x - \xi{Pi - &) + \x 2 {3f3 2 2 - 2fcfo - ft) + . . . = 0. 
These hodograph equations provide us with the solutions of the system 

In terms of the variables u — — (/3i + /? 2 ),w = — P2) 2 one has the dToda 
system (2). Considering the higher times £2,23, ■•■ -one gets the whole dToda 
hierarchy. 

Similar to the Benney case the function Wt is the generating function for 
classical singularities for functions of two variables. Indeed, in the variables 
X = §(& + /3 2 ), Y = \{fo- f3 2 ) it is of the form 

W T = -x X - \x x Y 2 - l -x 2 XY 2 - \x z {AX 2 + Y 2 )Y 2 + ... (39) 

The third term here represents the parabolic umbilic singularity both for hy- 
perbolic and elliptic cases. 

The formulas for the dToda hierarchy presented here coincide with those 
given in the paper [21] after the identification 

V T (x, X) = -2TX + AV£(t, A). (40) 

i.e. xq = —2T, x n = nt n , n = 1,2, 3, .... 

It is obvious that the descriptions of the regular and singular sectors of the 
dToda hierarchy completely coincide with those of 1-layer Benney hierarchy. 



10 



6 6. Interrelations between the Euler-Poisson- 
Darboux equations with different indices and 
those for function W and densities of integrals 
of motion. 



l-layer Benney hierarchy and dToda hierarchy are two examples of hydrody- 
namical type systems for which functions W obey the Euler-Poisson-Darboux 
equations 

d 2 



L F W F 



(-- 


--)} 




dy)\ 



W £ = 0, (41) 



-dxdy x — y 

with different indexes e. Such linear equations are well studied ( see e.g. 
[8]). The operators L £ have a number of remarkable properties. One of them ( 
probably missed before) is given by the identity 



Lfj, = L^ +1 L £ (42) 



for arbitrary indices e and fi . This identity implies, for instance, that for 
any solution W £ the function L^W £ with arbitrary \i obeys the Euler-Poisson- 
Darboux equation with index s + 1, more precisely L^W £ — e(e — jj,)W £+ \ . 
In particular, at e = — \ and u = one has LiL n = L\L_i. In terms of the 

operators L £ defined as L £ = (x — y)L £ the last relation takes the form 

d x d y L_i = Lid x d y . (43) 

This identity clearly demonstrates the duality between the Euler-Poisson-Darboux 
equations with indices \ and - \ and consequently between l-layer Benney and 
dToda hierarchies. 

Duality between the functions W and densities of integrals of motions is the 
another type of duality typical for the so-called e integrable hydrodynamical 
type systems. Indeed, due to the Tsarev's result [22] , a symmetry Wi of a 
scmi-Hamiltonian hydrodynamical system 

f = A, OT f, i = !,...,„, (44) 
i.e. a solution of the system 

97 = Wj(/ ^' t = 1 '-' n ( 45 ) 
which commutes with the system (44), are defined by the system 



dw k d\k 

dpi d/3 z 



Wi — wk Xi — Afe ' 



i ^ fc. (46) 
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Such Wi provide us with the solutions of the systems (44) via the hodograph 
equations 

fl t := -x + \i(j3)t + Wi =0,i= l,...,n. (47) 
For such system (44) densities P of integrals of motion obey the equations [22] 

dp % dp k \ t -\ k dp k ' T ' y ' 

Let us define e-systems as those ( for particular class of such systems see e.g. 
[23]) for which 

(49) 



Aj — Afe Aj — Afe Pi — pk 
For such systems densities of integrals obey Euler-Poisson-Darboux equations 

d 2 P e dP e dP 



dpidp k p l ~p k dp l p l ~p k dp k 

At the same time the equations for Wi become 



i ^ k. (50) 



dw k _ Wi - w k 

dP, " £ p^p k ' (51) 

Symmetry of these equations with respect to the transposition of indices i and 
k implies that = f§£- Hence 

dW 

W i = -^T, i = l,...,n, 

dp t 

for a certain function W. Thus, equations (51) arc the Euler-Poisson-Darboux 
equations of the type E(—e, — e) for the function W 

d 2 W / e dW £ dW 



dp t dp k 



f £ OW _£ CW\ 

\fr-Pkdp* Pi-p k dp k )' ^ [bZ) 



The fact that the generating function for symmetries of the Whitham equations 
and some other integrable hydrodynamical systems obey the Euler-Poisson- 
Darboux equations has been observed earlier in the papers [9, 13, 23]. Also 
the duality between the Euler-Poisson-Darboux equations for the densities of 
integrals of motions and generating functions of symmetries has been noted be- 
fore too. However the demonstration presented above seems to be different from 
those discussed earlier. 

In addition one can note that equations (49) imply that for e-systems also 
A, = with some function g. As the result the hodograph equations (47) for 
the e-systems take the form 

dg dW dW 

dPi dPi dPi 
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where W = —x(0i + P2) +gt + W. Thus, hodograph equations for the integrable 
hydrodynamical type systems are nothing but the equations defining the critical 
points of the function W. It seems that this fact has been missing in the 
previous publications. Moreover, due to the equations (49) the function g also 
obeys the E(—e, — e) Euler-Poisson-Darboux equation and , hence, the function 
W does the same. Note that particular class of e-systems for which \ are linear 
functions of /3j has been discussed in [23]. 

So, for integrable hydrodynamical systems of the e type, the densities of 
integrals and the functions W ( as well as the functions W generating sym- 
metries) play a dual role obeying the Euler-Poisson-Darboux equations with 
opposite sign of the index e . This property resembles a lot the well-known 
duality between the generating functions of integrals of motion and symmetries 
for the dispersionful integrable equations. 
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